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Abstract 

In this paper, we derive the exact rate of convergence of some approximation schemes 
associated to scalar stochastic differential equations driven by a fractional Brownian 
motion with Hurst index H. We consider two cases. If H > 1/2, the exact rate of 
convergence of the Euler scheme is determined. We show that the error of the Euler 
scheme converges almost surely to a random variable, which in particular depends on 
the Malliavin derivative of the solution. This result extends those contained in [T7] 
and [18]. When 1/6 < H < 1/2, the exact rate of convergence of the Crank-Nicholson 
scheme is determined for a particular equation. Here we show convergence in law of 
the error to a random variable, which depends on the solution of the equation and an 
independent Gaussian random variable. 
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1 Introduction 



Let B = (Bt, t £ [0, 1]) be a fractional Brownian motion (in short: fBm) with Hurst 
parameter H E (0,1), i.e., B is a continuous centered Gaussian process with covariance 
function 

R H (s,t) = ^(s 2H + t 2H - \t - s\ 2H ), s,t€ [0,1]. 

For H = 1/2, B is a standard Brownian motion, while for H ^ 1/2, it is neither a semi- 
martingale nor a Markov process. Moreover, it holds 

(ElBt-Bsl 2 ) 1 / 2 = \t-s\ H , a,te [0,1], 

and almost all sample paths of B are Holder continuous of any order a £ (0, H). 

In this paper, we are interested in the pathwise approximation of the equation 

X t = x + [ a(X s )dB s + f b(X s )ds, t e [0, 1], (1) 
Jo Jo 

with a deterministic initial value xq € R. Here, a and b satisfy some standard smoothness 
assumptions and the integral equation (JTJ) is understood in the sense of Russo-Vallois. Let 
us recall briefly the significant points of this theory. 

Definition 1 (following [21]) Let Z = (Zt)t£{o,i] be a stochastic process with continuous 
paths. 

• A family of processes (Ht^)tE[o,i] ^ s sa ^ t° converge to the process (Ht)te[o.i] * n ^ e UC P 
sense, i/sup te [ 01 ] \h[ 6 ^ — H t \ goes to in probability, as e — ► 0. 



The (Russo-Vallois) forward integral L Z s d B s is defined by 



lim-ucp e 1 / Z t (B t+£ - B t )dt, 
Jo 



(2) 



provided the limit exists. 
• The (Russo-Vallois) symmetric integral J Z s d°B s is defined by 

lim-ucp (26)" 1 I {Z t+£ + Z t ){B t+e - B t )dt, (3) 
Jo 

provided the limit exists. 

Now we state the exact meaning of equation ([T]) and give conditions for the existence and 
uniqueness of its solution. We consider two cases, according to the value of H: 
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Case H > 1/2. 

Here the integral with respect to B is defined by the forward integral ([2]). 

Proposition 1 If a 6 C? and if b satisfies a global Lipschitz condition, then the 
equation 

X t =x + [ a{X s )d-B s + [ b(X s ) ds, t £ [0, 1] (4) 
Jo Jo 

admits a unique solution X in the set of processes whose paths are Holder continuous 
of order a > 1 — H . Moreover, we have a Doss-Sussmann type J?[ \29$ representation: 

X t = <P(A t ,B t ), 4 6 [0,1], (5) 

where 4> and A are given respectively by 

- — (xi,x 2 ) = a((f>(xi,x 2 )), 0(xi,O)=xi, xt,x 2 eR (6) 
ox 2 

and 

4 = exp(- a\(j>{A t ,s))d^Jb{<t>(A u B t )), A = x , te[0,l}. (7) 

Proof. If X and Y are two real processes whose paths are a.s. Holder continuous 
of index a > and j3 > with a + j3 > 1, then J^Y s d~X s coincides with the 
Young integral J^YgdXg (see [28], Proposition 2.12). Consequently, Proposition [1] is 
a consequence of, e.g., [H] or [23] . 

□ 



Case 1/6 < H < 1/2. 

When H < 1/2, in particular the forward integral Jq B s d~ B s does not exist. Thus, 
in this case, the use of the symmetric integral ([3]) is more adequate. Here we consider 
only the case 6 = 0: for the general case see |19j . |21j and Remark [1] 

Proposition 2 // H > 1/6 and if a £ C 5 (M) satisfies a global Lipschitz condition, 
then the equation 

X t = x + [ a(X s )d°B s , t G [0, 1] (8) 
Jo 

admits a unique solution X in the set of processes of the form Xt = f(Bt) with 
f £ C 5 (R). The solution is given by X t = 4>{xo,B t ), t £ [0,1], where (f> is defined by 



Proof. See p2], Theorem 2.10. 



□ 



Remark 1 In [21], Nourdin and Simon developed recently a new concept, namely the 
Newton-Cotes integral corrected by Levy areas, in order to study equation (pQ) for any 
H G (0, 1). It allows to use a fixed point theorem to obtain existence and uniqueness in the 
set of processes whose paths are Holder continuous of index a G (0, 1) and not only in the 
more restrictive - and a little arbitrary - set of processes of the form Xt = f(Bt,At) with 
/ : 1 x [0, 1] -» 1 regular enough and A a process with (^-trajectories, as shown in [19] . 



Approximation schemes for stochastic differential equations of the type (pQ) are stud- 
ied only in few articles, see, e.g., [17J and the references therein. In [18], the second-named 
author considers the approximation of autonomous differential equations driven by Holder 
continuous functions (of any fractal index < a < 1). He determines upper bounds for the 
order of convergence of the Euler scheme and a Milshtein-type scheme, see also [30], and 
applies then his results to the case of the fBm. In [T7], the first-named author studies the 
following equation with additive fractional noise 

X t = x + I (j{s)dB s + [ b(s,X s )ds, t€[0,l] (9) 
Jo Jo 

under the hypothesis H > 1/2. For a mean-square-L 2 -error criterion, he derives by means 
of the Malliavin calculus the exact rate of convergence of the Euler scheme, also for non- 
equidistant discretizations. Moreover, the optimal approximation of equation @ is also 
studied in [T7] . 

In this paper, we are interested in the exact rate of convergence of the Euler scheme 
associated to @ and of the Crank-Nicholson schemes associated to ([8|). Thus here, com- 
pared to [17] . we study the non-additive case. We obtain two types of results (see Section 
3 for precise statements): 

(i) If H > 1/2 and under standard assumptions on a and b, then the classical Euler scheme 
X with step-size 1/n for equation defined by 

3 = X ° — n —n — « 1 (10) 

X (k+l)/n = X k/n + v( X k/n)AB k / n + b(X k/n )-, k € {0, . . . , H - 1}, 

and X t = Xt nt y n for t G [0, 1] verifies 

n 2H-l ^ _ Xi ] a* _ 1 I a ^ Xs)DsXlds _ (n) 

1 Jo 

Here D s Xt, s, t G [0, 1] denotes the Malliavin derivative at time s of Xt with respect to 
the fBm B. This result is somewhat surprising because it does not have an analogue 
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in the case of the standard Brownian motion. Indeed, in this framework, or more 
generally when SDEs driven by semimartingales are considered, it is generally shown 
that X 1 converges a.s. to X% and then that the correctly renormalized difference 
converges in law, see, e.g, [13] and Remark O point 3. For the approximation of Ito- 
SDEs with respect to mean square error criterions, see, e.g., [12j or [TBj and Remark 
[21 point 3. 

Moreover, if we consider the global error on the interval [0, 1] of the Euler scheme, we 
obtain 

a.s. 1 



n 2H - x ipr 



X 



z te[o,i] 



a\X s )D s X t ds 



(ii) Assume that 1/6 < H < 1/2 and let us consider the Crank-Nicholson scheme X n with 
step-size 1/re associated to (JS]): 



A,', 



Xq 



yn 

^{k+l)/n 



X k/n + 2 



2 (°"(^fe/n) + a (X(k+l)/S) ( B (k+l)/n ~ B k/n ), 

k G {0, . . . , n — 1}, 



(12) 



and A™ = XP t y n for i G [0, 1]. Here we obtain the following rates of convergence: 
• (Exact rate) If the diffusion coefficient a G C 1 satisfies 

<t(x) 2 = qx 2 + /?£ + 7 with some a, /?, 7 G M, 

we have 



n 



3ff-l/2 



A" - Ai 



(13) 



with G ~ N(0, 1) independent of Ai and a\ given by (f22|) . We prove also an 
equivalent of (fT3|) . at the global level: 



n 3H ~ 1/2 sup 

fee{o,...,n} 



V n 



A 



fc/n 



— sup |o-(A() W t | 



12 



te[o,i] 



with W a standard Brownian motion independent of A. 

Compared to the above result for the Euler scheme, the convergence to a mixing 
law obtained here is classical in the semimartingale framework. In the fBm 
framework, such a phenomenon was already obtained in three recent papers for 
1/2 < H < 3/4: in [5], the authors study the asymptotic behavior of the power 
variation of processes of the form u s dB s , while, in [13] and [15] the asymptotic 
behavior of J Q * /(A^)G(A™ n^ _1 )ds is studied, where A™ denotes the broken-line 
approximation with stepsize 1/n of the solution A of (JSj) and A™ its derivative. 
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• (Upper bound) If 1/3 < H < 1/2 and a £ is bounded, we have 

for any a < 3H - 1/2, n a - ^ 0. (14) 

Note that (i) covers in particular the case of a linear diffusion coefficient, while in (ii) 
we consider smooth diffusion coefficients, which are bounded and therefore nonlin- 
ear. The exact rate of convergence in the general case (ii) remains an open problem 
although it seems that it is again 3H — 1/2, as in (i). 

The paper is organized as follows. In Section 2, we recall a few facts about the Malliavin 
calculus with respect to the fBm B. Section 3 contains the results concerning the exact 
rates of convergence for the Euler and the Crank-Nicholson schemes associated to and 
([8]) respectively. The proofs of the results for the Euler scheme are postponed to Section 4. 



2 Recalls of Malliavin calculus with respect to a fBm 

Let us give a few facts about the Gaussian structure of fBm and its Malliavin derivative 
process, following Section 3.1 in [23J and Chapter 1.2 in [22]. Let £ be the set of step- 
functions on [0, 1]. Consider the Hilbert space TL defined as the closure of £ with respect to 
the scalar product 

{\o,t]A[o,8]) H = RH(t,s), s,i€ [0,1]. 

More precisely, if we set 

K H (t, s) = T(H + l/2)-\t - sf- 1/2 F (H - 1/2, 1/2 -H;H + 1/2, 1 - t/s) , 

where F denotes the standard hypergeometric function, and if we define the linear operator 
K* H from £ to L 2 {[0,1}) by 

(K* H <p)(s) = K H (T,s)<f(s) + J ( v (r)-<p(s))-^(r,s)dr, <peH, se[0,l], 

then TL is isometric to L 2 ([0, 1]) due to the equality 

(<P,P) H = [ T (K* HV )(s)(K* H p)(s)ds, V ,peH. (15) 
J o 

The process B is a centered Gaussian process with covariance function Rh, hence its asso- 
ciated Gaussian space is isometric to TL through the mapping lr t i i— > B t . 

Let / : W 1 — > R be a smooth function with compact support and consider the 
random variable F = f{Bt 1 , . . . , Bt n ) (we then say that F is a smooth random variable). 
The derivative process of F is the element of L 2 (Q,TC) defined by 

DsF = ^2g^(B tl ,...,B tn )l [0:ti] (s), a €[0,1]. 

i=l 1 
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In particular D s Bt = lro i t](s). As usual, B 1 ' 2 is the closure of the set of smooth random 
variables with respect to the norm 

||F|| 2 j2 = E[|F| 2 ] +E 

The divergence operator 5 is the adjoint of the derivative operator. If a random variable 
u G L 2 (Q.,Tl) belongs to the domain of the divergence operator, then S(u) is defined by the 
duality relationship 

E(F5(u)) = E(DF,u) n , 

for every F G B 1 ' 2 . 

Finally, let us recall the following result proved in [20] : 

Proposition 3 Let H > 1/2, a G C 2 and b G C\. Then we have for the unique solution 
X = (X t , t G [0, 1]) of equation that X t G B 1 ' 2 for any t G (0, 1] and 

D s X t = a(X s )exp ^ jf b'(X u )du + o-'(X M )d _ B u ^ , < s < t < 1, i > 0. (16) 




3 Exact rates of convergence 

In the sequel, we will assume that [0, 1] is partitioned by {0 = to < h < . . . < t n = 1} 
with tk = k/n, < k < n. Rates of convergence will thus be given relative to this partition 
scheme. For simplicity, we write AB k / n instead of Bn.+\y n — B k / n . For the usage of 
non-equidistant discretization in the simulation of fBm, see Remark [3j point 2. 

3.1 Euler scheme 

In this section, we assume that H > 1/2 and we consider equation i.e., the integral 
with respect to B is defined by ([2]). The Euler scheme X for equation is defined by 
(jlOp . The following theorem shows that the exact rate of convergence of the Euler scheme 
is in general n~ 2H+1 for the error at the single point t = 1. 

Theorem 1 Let b G C 2 and a G C^. Then, as n — ► oo, we have 

n ™-i [X« _ Xl ]^-l f 1 a'(X s )D s X 1 ds. (17) 

1 Jo 

Proof. We postpone it to Section 4. n 
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Remark 2 1. The asymptotic constant of the error does not vanish, e.g., for the linear 
equation with constant coefficients, i.e., 



X t = xo + -r [ X s (TB s + j3 f X s ds, te[0,l] 
Jo Jo 



(18) 



with 7, [5 G M, 7 / and x ^ 0. 



2. The appearance of the Malliavin derivative in the asymptotic constant seems to be 
due to the fact that (D t Xi) te ^ 01 ^ measures the functional dependency of X% on the 
driving fBm, see |25| . 



3. For the Ito-SDE, i.e., the case H = 1/2, it is shown in [13] that 



n 



1/2 



Tl - x\ t& 



as n — > oo, with a Brownian motion W, which is independent of the Brownian motion 
B, and 



Y s = exp U S b'(X l f) - \aa'{X™) du + jT a'{Xf>) dB u j , s G [0, 1]. 
In [3] the analogous assertion for the mean-square error is established, i.e., 

2 



nE 



Tl - X\ t& 



1 



E 



Fi / aa'{Xf 6 )Y~ l dW t 



as n — > cxd. 



4. Assume that 6 G C^, a G and that additionally 6 and cr are bounded with 
inf zg K |<r(x)| > 0. Under these stronger assumptions, which are a priori only of 
technical nature, we can show - by applying the same techniques - that Theorem [T] is 
also valid with respect to the mean square error, i.e., 



n 2H ~ l (eIA^-XJ^ 1 " 



E 



f a'{X s )D s X 1 ds 
Jo 



2\ V2 



as n — ► oo. Note that under the above assumptions on b and a equation simplifies 
to 

b(0(A t ,B t )) 

, a = x , te[0,i\, 



A' 



1 a{4>(A u B t )) 

which allows us to control the integrability of the remainder terms in the error expan- 
sions made in the Proof of Theorem [TJ 
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5. For equation (|18p and when H < 3/4, it is possible to go further and to obtain a 
convergence in law for the third term in the asymptotic development of X l , see also 
Theorem EJ We have, as in +00, 



x; ^ x 1 



n 2H-l/2 



7! v. «1-2H 



X'[ -X x + \X x n 



a.s. 



n 2H-l ^ _ Xi 

r 7 2 



Xi (19) 



with G a centered Gaussian random variable. For H = 3/4 the last convergence is 
again valid if one replaces n 2H ~ 1 / 2 by n 2H ~ 1 ^ 2 (logn) -1 / 2 . Indeed, we have 



n-1 



X? = x exp Lb 1+ P- i 7 2 ^(A£ fc/n ) 2 + R n 



k=0 



with n^-^l^l^o as n — > 00. Thus it holds, as n — > +00, 



2 2 n— 1 

fc=0 

Hence Theorem 3 (or Theorem 6 for the case H = 3/4) in [5] allows us to obtain the 
convergence in law in (|19p . When H > 3/4, it seems to be hard to derive a result in 
law since, in this case, arguments used in the proof of Theorem 3 in [5] are not valid 
anymore. Indeed, in this case, we do not work in a Gaussian framework, see, e.g., 
Theorem 8 in [5]. 

To overcome this problem one can modify the Euler scheme for the linear equation 
such that the second order quadratic variation of B appears in the error expansion. 
See, e.g., [H] for a similar strategy in the case of weighted p-variations of fractional 
diffusions. The second order quadratic variation of fractional Brownian motion is 
given by 

n-1 

\2 



k=l 

It is well known, compare for example [T], that 

n 2H ' l V 2 (B) ^ A-2 2H 
n 2H - l / 2 V 2 {B)-n 1 / 2 {A-2 2H ) -±> G, 

as n — ► 00, where G is a centered Gaussian random variable with known variance 
c 2 H > 0. Moreover, by an obvious modification of Proposition d] we also have 

B x ,n 2H - x l 2 V 2 (B) - n^ 2 (4 - 2 2H )) (Bi, G) 



9 



as n — > oo, with G independent of B\. 

For the following approximation scheme for the linear equation (fl"8 



X ™ = x 

X ( k+ l)/n = X k/n + l X k/n^k/nB + $Xj> /n A U/n B A {k _ 1)/n B + (3X^, 



n I 

k/nn' 

k £ {0, . . . ,n - 1}, 



we get by straightforward calculations 

X? = X x exp ("-I 7 2 v; 2 (i?) - - f3 1 B 1 - + i? n 
\ 4 An n 

with n mm { 2 ^ -1 / 2,1 }|-R n | -^>0 as n — ► oo. Hence we obtain 

„2 



-(4-2 2W )Xi 



as n — ► oo. Furthermore, we get the following error expansions according to the 
different values of H . 



(i) Case 1/2 < H < 3/4: 
2H-1/2 



n 



X 1 n -Xi + ^(4-2 2ii )X 1 n 1 



>2i?\ v 1-2H 



(ii) Case # = 3/4: 



Xf - Xi + -L(4 - 2 3 / 2 )X 1 n~ 1 



/2 



£ 7 2 /? 2 



(iii) Case 3/4 < H < 1: 



n 



X 1 ri -Xi + ^(4-2 2ii )Xm 



l-2# 



/? 2 



Thus for this scheme we get - according to the values of H — different error expansions 
due to the drift part of the equation. 



For the global error on the interval [0, 1], we obtain the following result. 
Theorem 2 



n^pT-XU^ 1 - sup 

1 te[o,i] 



<j'{X s )D s X t ds 



(20) 



Proof. We postpone it to Section 4. 
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□ 



Hence the Euler scheme obtains the same exact rate of convergence for the global error on 
the interval [0, 1] as for the error at the single point t = 1. Moreover we have a.s. 



if and only if a.s. 



sup 

te[o,i] 







f a'(X s )D s X 1 ds 
Jo 

X t G (ff<r / )~ 1 ({0}) for all t G [0, 1]. 



Remark 3 1. If b = 0, Theorem [T] and [2] are again valid under the weaker assumption 
that a G C\ . Since in this case 

D s X t = a{X s )exp ( I a'{X u )d-B u ) = a{X t ) s G [0, t], 



which is an obvious consequence of the change of variable formula for fBm, we have 
here 

^ [Tl -X,}^- \a{X x ) I'' a'(X s )ds 



n 



and 



respectively. 



n 



2H-1 



\X -X\\ sup 

Moo o r r , 

Z *6[0,1] 



Jo 



2. For H ^ 1/2 the increments of fractional Brownian motion are correlated. There- 
fore the exact simulation of B(ti), . . . ,B(t n ) is in general computationally very ex- 
pensive. The Cholesky decomposition method, which is to our best knowledge the 
only known exact method for the non-equidistant simulation of fractional Brownian, 
requires 0(n 3 ) operations. Moreover the covariance matrix, which has to be decom- 
posed, is ill-conditioned. If the discretization is equidistant, i.e., = i/n, i = 1, ... ,n, 
the computational cost can be lowered considerably, making use of the stationarity of 
the increments of fractional Brownian motion. For example, the Davies-Harte algo- 
rithm for the equidistant simulation of fractional Brownian motion has computational 
cost 0(n log(n)), see, e.g., [6]. For a comprehensive survey of simulation methods for 
fractional Brownian motion, see, e.g., [4]. 

3. In the Skorohod setting, it is in general difficult to write an Euler type scheme asso- 
ciated to equation ([T]), even if H > 1/2 and 6 = 0. Indeed, in this case, by using the 
integration by parts rule 5(Fu) = FS(u) — (DF, u)-^ for the Skorohod integral and by 
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approximating X (fe+1)/n by X k/n + J k ( /n >/n a(X k/n )5B s (as in the case H = 1/2), one 
obtains 

X {k+l)/n = X k/n + a i X k/n) { B (k+l)/n ~ B k / n ) - a'(X k / n ) {DX k / n , l[k/n,(k+l)/n])n- 

The problem is that the Malliavin derivative DX k / n appears, which is difficult to 
compute directly. Moreover, the error analysis of such an approximation seems also 
to be very difficult, because the L 2 -norm of the Skorohod integral involves the first 
Malliavin derivative of the integrand. Thus, for analyzing such an approximation 
scheme, we need also to control the difference between the Malliavin derivative of 
the solution and the Malliavin derivative of the approximation. But this involves the 
second Malliavin derivative etc. and we cannot have closable formulas. It is one of 
the reasons for which we preferred here to work within the Russo-Vallois framework, 
instead of the Skorohod one. Another reason is that the Russo-Vallois framework is, 
from our point of view, simpler in the one-dimensional case than the Skorohod one, 
as it is shown in |19j . 



3.2 Crank-Nicholson scheme 

In this section, we assume that 1/6 < H < 1/2 and we consider equation ([8]). Let X n 
be the Crank-Nicholson scheme defined by (|12p . which is the canonical scheme associated 
to equation ([8]), since the integral with respect to the driving fBm B is defined by the 
symmetric integral. It is an implicit scheme, but it is nevertheless well-defined, since for 
n sufficiently large, ihi- ^AB k / n <r(x) is invertible. Although (1121) seems to be rather 
close to pop with 6 = 0, the situation is in fact here significantly more difficult. That is 
why we study the rate of convergence for the Crank-Nicholson scheme only in the following 
particular cases: 

• Case 1: 1/6 < H < 1/2 and a € C 1 satisfies a(x) 2 = ax 2 + [3x+^f for some a, j3, 7 E R, 

• Case 2: 1/3 < H < 1/2 and a G C fe °° bounded. 



3.2.1 Case 1 

Compared to TheoremHJ we have here a convergence in law. Moreover, the limit of the error 
is expressed as a mixed law between B\ and an independent standard Gaussian random 
variable G, see also Remark [21 point 4. 



Theorem 3 Assume that 1/6 < H < 1/2 and a £ C 1 satisfies a(x) 2 = ax 2 + (3x + 7 for 

some a, /3, 7 G M. Then, as n — > 00, we have 



n 



3H-1/2 



x? - x l 



^an^aiX^G. 



(21) 



12 



Here G ~ N(0, 1) is independent of Xi and 

oo 

a 2 H = 4/3 + 1/3 9 ( l f> where 2e ( £ ) 



+ 1 



,2H 



(t-1 



} 2H _ 2£ 2H 



(22) 



In fact, we have also a result at a functional level: 



Theorem 4 Assume that 1/6 < H < 1/2 and a € C 1 satisfies a{x) 2 = ax 2 + (3x + 7 /or 

some a,/?, 7 £ R. Then, as n — > 00, we have 



n 3H ~ 1/2 sup 

ke{0,...,n} 



X k/n ~ X k/n 



— >0 ~H— sup \a(X t )Wt\ 



12 



te[o,i] 



(23) 



i7ere W is a standard Brownian motion independent of X and cth > is once again given 
by (El). 



Remark 4 1. In [19|, equation ([8]) is also studied and it is shown that X™ converges in 
probability if and only if H > 1/6 and that, in this case, the limit is X\, the solution of 
equation (|8|) at t = 1. Of course, this fact is also an obvious consequence of Theorem 

m 

2. Let us show how the constant in (f22l) appears. Set Y n = n~ 1 / 2 Ylk=o( nH ^-Bfc/n) 3 . We 

claim that Var[y n ] -> a 2 H as n -> +00. Indeed, using (n H 5 t / n ) tg [ 0jOO ) = (B t )te[o,oo), 
we have 

n-l 

Var[F n ] = n- 1 £ E [(B fc+1 - B k ) 3 (B e+1 - B e ) 3 ] . 
k,e=o 

Since x 3 = H%(x) + 3i?i(x), where H\(x) = x and H^(x) = x 3 — 3x denote the 
first and third Hermite polynomial, we can write, by using the well-known identity 
E[Hi(X)Hj(Y)] = for i ^ j and E[Hi(X)Hi(Y)] = E[XYf/i\ for a centered Gaus- 
sian vector (X,Y) with Var[X] = VarfK] = 1 (see, e.g., [22], Lemma 1.1.1): 

n— 1 n— 1 

Var[y n ] = (6n)- x £ - /c) 3 + 9^ £ - k), 
k,e=o k,e=o 

with 



29(£ - k) = 2E [(B fc+a - - B t )\ = \£ - k + 1\ 2H + \£-k- 1\ 2H - 2\l - k\ 

On the one hand, let us remark that n~ l J2k~e=o = n~ l E (j2k=o ^ 

n 2H ~ x — > Oasn-> +00, if H < 1/2. On the other hand, we can write 



n— 1 Ti—l 7i—l n—1 n— In— fc— 1 

e(£-k) 3 = Y,0(0) 3 + 2Y^ £ ^-fc) 3 = 8n + 2^ £ 

fe,<2=0 fc=0 fe=0^=fc+l fc=0 £=1 
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Consequently, since 6(1) < for H < 1/2, we deduce by using Cesaro's theorem that 
Var[F„] — > a\ given by (f22j) . as n — > +oo. 

3. In [10] , in particular the approximation of fBm by its piecewise linear interpolation 

Bt = B k/n + (nt - [nt)){B {k+l)/n - B k/n ), t G [k/n, (k + l)/n] 

is studied. It is shown that the correct renormalization of \\B — B n \\ OQ converges to 
the Gumbel distribution, i.e., 

P [\\B n - BWoo < a n (v n + x/v n )^ — > exp(-exp(-x)) 

as n — > for i6R, where a n ~ cnn~ H with ch > and v n is in terms of log(ra). Since 
3H — 1/2 > H for .ff > 1/4, the analogue of Theorem 2 in the setting of Theorem 4, 
i.e., 

£ a 



^H-l^^^J^^ gup 1^)^ 

12 te[o,i] 

as n — > oo, can not hold without further restriction of the Hurst parameter. 



For the proof of Theorem [3l we need the following Lemma. 
Lemma 1 i) We have, for H < 1/2, 



(V^- 1/2 ^(AB fc/n ) 3 J -^(B 1 ,G), aS n^+oo, (24) 
V fc=0 / 

where G is a centered Gaussian random variable with variance a\ given by (f£2|). 
independent of B\ . 



ii) We have 



n-1 

n 

k=0 

where G' is a centered Gaussian random variable. 
Hi) We have 

n-1 

6ff-lVV a d . \6 Prob 



^^(Ai?fc/ n ) 6 15, as n — > +oo. 



fc=o 



Proof of Lemma [TJ The second and the third point are classical: we refer to [2]. Thus 
we have only to prove the first point. Let us denote by H\{x) = x and H$(x) = x 3 — 3x the 
first and third Hermite polynomial. Since H < 1/2, we have 



^J2 Hl (n H AB k/n ) = n H -^B x ^ 



n 

k=0 
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and we deduce that the convergence in law (|24|) will hold if and only if 

(B 1 ,n- 1 ^H 3 (n H AB k/n ) \ -^(B 1} G), (25) 
V fc=o / 

as n — > +00. In [2], it is shown that n _1//2 Ylk=o Hz(n H ABf,/ n ) — — > G. See also the second 
point of Remark |U Then, the proof (|25p is finished by the following Proposition, which is 
an obvious consequence of the main result contained in [26]. □ 



Proposition 4 Let (F^,F^) be a random vector such that, for every n, Ff (i = 1,3) is 

c c 

in the i-th Wiener chaos associated to the fBm B. If F™ — > G\ and F£ — >G% with G{ 
(i = 1,3/ some Gaussian variables, then C(F™ , F£) — > £(Gi) <8> £-{Gz)- 



Proof of Theorem [3j In [19] . the second-named author proved 

(re— 1 2 n ~ 1 /re— 1 \ \ 

B X + ^ ^(A^/J 3 + ^ £(Ai? fe/n ) 5 + O ^(A5 fc/n ) 6 po ■ (26) 
fe=0 fc=0 \fe=0 / / 

Note that X\ = 4>{xq, B\) and Jjf-^i, x s ) = a{4>(xi,X2))- Consequently by a Taylor expan- 
sion Xf — X\ equals 

re— 1 2 n— 1 / re— 1 /re— 1 \ 2 

^a(X 1 )^(A J B fc/n ) 3 +^a(X 1 )^(A J B fc/n ) 5 +0 1^1 x ^(A^/J 6 A ^(Ai? fc/n ) 3 

A;=0 ' fc=0 \ fc=0 \fc=0 / 

By the second point of Lemma [Tj we deduce that n 3 ^- 1 / 2 Efc=o( A5 fc/n) 5 — ► in law, 
hence in probability. By the third point of Lemma [Tj and since H > 1/6, we have that 
^ 1/2 S=o( A ^/«) 6 = n l / 2 ~ m x n 6 "- 1 ££-(!( A5 fc/n ) 6 — ► in probability. By the 

first point of Lemma HJ and since H > 1/6, we have that n 3 ^" 1 / 2 ( X^A-=o( A -^fc/r 



2 



n l/2-3H f n 3H~l/2 S^?~h 



2 



Ylk=o(^--Bk/n) ) — > in probability. Then, using again the first point 
of Lemma [T]and Slutsky's Lemma, we obtain (|21j) . U 

Proof of Theorem [4[ Exactly as in the proof of Theorem 3 in [5] , we can prove 

(B u n ZH - l ' 2 Vi l (B) t ) -±> <r H W t ), as n - +00, 

in the space £>([0, l]) 2 equipped with the Skorohod topology. Here, W is a standard Brow- 
nian motion independent of B and V£{B) is defined by V^{B) t = ^[=o _1 (A^/ n ) 3 . To 
obtain Theorem [H it suffices then to adapt the proof of Theorem [3] below. r-. 
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3.2.2 Case 2 

Now assume that a G C£° is bounded and that 1/3 < H < 1/2. 

In the sequel, we will need some fine properties concerning the m-order variation of 
B on the interval [0, 1]. Let us state them in the following proposition: 

Proposition 5 Let h £ C^. 

1. If m £ N is even then, for any H £ (0, 1): 

n-l 



n mH- 



1 ^K B k/n) ( B (k+l)/n -B k/n ) m ^ 2^/2^/2)! /o h( - Bs " ,dS ^ asn ^+ co - 



k=0 

2. IfmeN \ {1} is odd i/ien, /or any i? G (1/4, 1/2) and a < ml - 1/2: 



n-l 



h( B k/n) + 2 h'(B k / n )(B( k+1 y n - B k / n ) 



(-S(fc+l)/n - Bk/n) m 0, (27) 



fc=0 

as n — ► +oo. 
Proof. 

1. When /i = 1, it is a classical result: we refer to [2] for instance. To obtain the general 
case, it suffices to adapt the methodology developed in step 5 of [8j, p. 8. or in the 
proof of Theorem 1 in [5] . 

2. Using the same linear regression as in the proof of Theorem 4.1 in [9], we can prove 
that, when H < 1/2: 

n-l 

n a y,[fc(ff(fc+l)/ n ) + H B k/n)\ (B(k+l)/n ~ B k / n ) m 
k=0 

converges in probability to 0, for any a < mH — 1/2. Convergence (|27p can then be 
obtained using a Taylor expansion and the fact that H > 1/4. 

Details are left to the reader. _ 

Using the above proposition we can show the following result. 
Theorem 5 Assume that H £ (1/3, 1/2) and a £ is bounded. Then we have: 

For any a < 3H - 1/2, n a \x[ n) - Zi] ^ (28) 

as n — > oo. 
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Proof of Theorem [5], In the following, denote 

A J z k / n = (Z( k+1 y n - z k / n y 

for j, n G N, k G {0, . . . , n — 1} and a process Z = (Z t )te[o t i]- When j = 1, we prefer the 
notation AZ k / n instead of A Z k i n for simplicity. Denote also, for p £ N, 

A%B)= max \A p B k/n \. 

fe=0,...,n— 1 



Moreover, it is recall that 4> given by © verifies the semigroup property: 

\/x,y,z£R, 4>(<j)(x,y),z) = $(x;y + z) 

and we have 

X t = <P(x,B t ), i€[0,l]. 
Simple but tedious computations (see, for instance, [19]) allow us to obtain: 

ltD/n = * (%%, A% + /(^i)A 3 % + g(X^ n )A*B k/n + 0(A\B))) , 



(29) 



with / = ( a „) and g 



every £ 6 {1, . . . , n}: 



12 



-. We deduce from the semi-group property (f29j) that for 



W =4f U Bljn + /(iW)A 3 B fc/n + £ 5 (lg)A 4 B fe/n + 0(nA 5 (B)) 



i-i 



e-i 



(30) 



fc=0 



fc=0 



In particular, we have 



sup 

te{o,...,n} 



sup 

£G{0,...,n} 



(n) 



and then (l30l) becomes 



(x,B i/ri/ 



0(nA 3 (B)) (31) 



^ = 4> U B i/n + Y, f(X ( k %)A 3 B k/n + Y, g(X k/n )A 4 B k/n + 0(n 2 A 7 (B)) . (32) 

V k=0 k=0 / 

Due to the assumptions on a we have that 

<f>(x,y 2 ) = <t>(x, yi ) + J2 Jijd^y^' Vl)( - V2 ~ Vl)j + ~ yi)m+1) - 
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Thus we get 

-H _ v . f 1 di<f> 



X k/n ~ X k/n + 2^ J JforfW Bk l n) 
fc-1 fc-1 

£ /(xg n )A 3 B fcl/n + ^ ? (l, l/n )A 4 B, l/n + 0(n 2 A 7 (B)) 

v fci=0 fei=0 
+0(n m+1 A 3 ( m+1 )( J B)). 

(i) Now assume for a moment that H > 5/12. By using (f33l) with m = 1 and ^^(a?i,X2) 
(t(0(xi,jc 2 )) we get 



(33) 



fc-i 

l fc/n 

fci=0 



X™ = X k/n + a(X k/n ) f(Xl /n )A 3 B kl/n + 0(n 2 A e (B)). 

and 



X™ = X k/n + a(X k/n ) Y, f(X kl/n )A 3 B kl/n + 0(n 2 A 6 (B)). (34) 

fci=o 

By inserting the previous equality in (|32[) with £ = n, we obtain 

(n— 1 n— 1 fc 

B l + f(X k/n )A 3 B k/n + ^ /V(X fe/n )A 3 i? fe/n ]T /(X fel/n )A 3 B fel/n 
fc=0 fc=0 fci=0 

n-1 \ 

+ 9(X k/n )A*B k/n + 0(n 3 A 9 (£)) . (35) 

k=0 J 

But we have, due to the second point of Theorem [5] and the fact that g = 2 af and 

X t = 0(x,B t ): 

(n-l n-1 \ 

^/(X fc/n )A 3 J B fc/n + ^ 5 (X fc/n )A 4 S fc/n ^5 0. 
fc=0 fc=0 / 

Moreover, since H > 5/12 and a < 3H — 1/2, we have 

n a+3 A 9 (5) ^ 0. 

On the other hand, since 

n— 1 A; 

2 ^ f' (T ( x k/n)A 3 B k / n y f(X kl / n )A 3 B kl / n 

k=0 fci=0 

(n-1 \ / n-1 \ n-1 

Y,f'v(X k/n )A 3 B k/n ) f(X kl/n )A 3 B kl/n \ - ^ / A(X fc/n )A 6 £ fc/n 

fc=0 / \fci=0 J k=0 

18 



we deduce, this time due to the first and the second point of Theorem [5l that 

n— 1 k 
n a J2f'°( X k/n)^B k/n £ f(X kl/n )A 3 B kl/n ^ 0. 
fc=0 fci=0 

Finally, we obtain ([25]) when H > 5/12. 

(ii) To prove the announced result, that is (j28l) for arbitrary H £ (1/3,1/2), it suffices to 
use ([33]) with the appropriate m for the considered H and then to proceed as in (i). The 
remaining details are left to the reader. 

□ 



4 Proof of Theorems [T] and [2] 

Throughout this section we assume that b E C^, a G Cf and H > 1/2. For g : [0, 1] — ► K 
and A G (0, 1) we will use the usual notations 

n n | ( .s IffO) -g(s)\ 

1151100= sup \g(t)\, \\g\\\= sup — n- — . 

te[o,i] s,te[o,i],s^t \ l ~ s \ 

Moreover positive constants, depending only on b, a, their derivatives, xq and H, will be 
denoted by c, regardless of their value. We will write A instead of 1/n. 



The following properties of the function </> are taken from Lemma 2.1 in [30J. 

Lemma 2 Let <f> given by |jS[). Then we have 

(a) (j)(x 1 ,x 2 ) = (p(4>(x 1 ,y),x 2 - y), x\ = <K<K^i, x 2 ), -x 2 ) 

(b) - — (x 1 ,x 2 ) = a(xi)- — (xi,x 2 ) 

OX 2 OX I 

(c) 0-2 fal) ^2 _X 2) - 2<7 ( Xl ) Q x Q x ( XU ~ X2 "> + ( aa ')( X l)fy^( X l> ~ x 2) 

(d) I = ^-{4>{xi,x 2 ) 1 -x 2 )^-{x 1 ,x 2 ) 

ax\ ox\ 

d(f> ( f X2 \ 

(e) — (a3i,s 2 ) = exp I y a ((f>(xi, s)) ds J 

for all x\,x 2 ,y G M. 



The following is well known and easy to prove. 
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Lemma 3 Let n G N and ai, b{ S R /or i = 1, . . . ,n. 

faj For Xj, jf = 0, . . . , n given by the recursion 

Xj + i = xjbj + aj, j = 0, . . . ,n — 1, 

mf/i = 0, we /icwe 



Xj = ^ «, ] [ Ofc, j = 1, ... 71. 

i=0 fc=i+l 

< + \a,j\, j = 0, . . . ,n - 1, 

luzf/i = and |&J > 1 /or all j = 0, . . . ,n — 1, then 



n— 1 n— 1 



max < |aj| TT |6 fc |. 

j=0,...,n •" 1 1 li. 1 1 



i=0 fc=l 



We will also require that the Euler approximation of the solution and the process (At)te[o,i] 
given by ([7|) can be uniformly bounded in terms of the driving fBm. 

Lemma 4 We have 

sup sup |X fc / n | < exp(c exp(c(||S|| 00 + ||5||w 2 ))) a.s., (36) 

neNfc=0,...,n 

Halloo < exp(c exp(c H-Blloo)) a.s. (37) 

Proof. We prove only the first assertion, following the proof of Lemma 2.4. in |30j . The 
second assertion can be obtained by a straightforward application of GronwalFs Lemma to 
equation © . By Lemma [2] (a) we have 

H^k/m -Bk/n) = HH^k/ni AB k / n ), — B ( k+1 )/ n ) . 

Using this, we obtain by the mean value theorem 

0(^(fc+l)/n> - B (k+l)/n)-4>( X k/ni _ ^fe/n) = 

with between / n and (j)(X k / n ,AB k / n ). Moreover 

n n n i n 

(t>{ X k/n^ B k/n) = 4>(X k/n ,0) + AB k/n -Q— (X k/n ,0) + -(AB k/n ) -Q^{X k/n ,Ck) 

with \( k \ < \AB k / n \. Thus we get 

<l>( x k/n, ^B k/n ) = Tl /n + a(X n k/n )AB k/n + ^(aa')(4>(X k/n , C, k ))(AB k/n f 



X 



(k+l)/n 



H X k/n,AB k / n ) 



■7T-(£k,-B {k+1)/n ) 
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by the definition of <p. So we finally obtain 



Since 



4>{ X (k+l)/ni —B(k+l)/n) ~ 4>{ X k/m ~Bk/n) 

b{l? k/n )A - ±(aa')(<p(X n k/n ,(k))(&B k/n ) 2 

F)A) ( f x 2 

-(xi,x 2 ) = exp I J a (4>(xi, s)) ds 



dx\ 



(6c> -B( k+1 y n ). 



by Lemma [2] (e) , we have 



d<j> 



(Cfc> -B(k+l)/n) 



<exp(c||B|| c 



Due to the assumptions, the drift and diffusion coefficients satisfy a linear growth condition, 
i.e., 

\b(x)\ <c{l + x), \a(x)\ < c(l + x) 

for x £ R. Hence we get 

\<t>(X(k+l)/n> -B(k+l/n))\ < \<P( X k/n, -B k /n)\ + c exp(c ||S|| 00 )(1 + |X£ /r J)A 

+c exp(c ||S|| 00 )(1 + \<j){Xl /w ( k )\)(AB k/n ) 2 . 

Since by Lemma [2] (a) 

X k/n = 4>{4>{ X k/w~Bk/n)^B k i n ), 

and (j)(0, 0) = 0, we have by Lemma [2] (b) and (e) 



\X k/n \ < c exp(c H^Hoo) [ \<l>{X k/n , -B k/n )\ + \\B\ 
and furthermore, since |Cfc| < \^B k j n \ 

\<t>(Xk/nXk)\ < c exp^HBlloo) (jXl /n \ + \\B\ 

< c exp(c H^Hoo) ( \cf>(X n k/n , -B k/n )\ + ||B[| C 

Together with 



\AB k/n \ < \\BWwA 1 ' 2 , 



this yields 



\4(X n {k+1)/n ,-B (k+1/n) )\ < \<f>(X n k/n ,-B k/n )\ [l + cexp(c\\B\U(l + \\B\\j /2 )A 

+cex P (c||S||oo)(l + ||S||J /2 )A. 
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Setting M = ||-B||oo + ||-6||x/2 ^ follows by LemmaH 

n n 

\(t>(X n k/n , -B k/n )\ < c exp(cM)A + c exp(cM)A) 

k=l j=l 

< exp(c exp(cM)). 

Thus with 

X k/n = 4>(4>{X k / n ,-B k / n ),B k / n ), 

we get the estimate 

\Xk/ n \ < ex P( c ll^lloo)|'/ , (^fc/n ) --Bfc/n)! + cexp(c H^Hoo) WB]^ < exp(c exp(cM)). 

□ 

Now we will state some Lemmas, which will be needed to determine the asymptotic constant 
of the error of the Euler scheme. The following Lemma [5] can be shown by straightforward 
calculations. 

Lemma 5 Denote 

f(x,y) = exp (- v\4>{y,s))ds\K4>{y,x)), x,y£R. 

Then we have f £ C 1,2 and in particular 

f x d<t> 

fy{x,y) = b'((p(y,x)) - f(x,y) J a (0(y,s)) — (y,s)ds, x,y £R. 
Lemma 6 We have a.s. 

ft 



< s < t < 1. 



exp QT b'{X u ) du + a'{X u )d-B u 

} -(A t ,B t )-^-(X s ,-B s )exp (J f y (B u ,A u )du 



dx 

Proof. By Lemma [2] (d) and (e) we have 



^-(A t ,B t )=exp^ ' (x'{<j){A t ,u))dv^, 
t\X 8 ,-B a ) = exp^-^ S u'{(j){A s ,u))duj 
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Using the notation 

g(x,y) = / a'(4>(y,u))du 
Jo 

we get by Lemma [5] 



dx 



'-{A t , B t )-^-(X s , -B s ) exp (J f y (B u , A u ) du 

= exp (g(B t ,A t ) - g(B s ,A s ) + J f y (B u ,A u )du) 
= exp (J b'{X u ) du^j exp (g(B t , A t ) - g{B Sl A s )) 

exp (- £ JJ U a '\<t>(A u ,T))^-(A u , r) dr f(B u , A u )d*?J . 



Since 

gx(x,y) = a'(4>(y,x)) 



and 



g y ( x ,y)=l v'\Hy> s ))-Q^(y> s ) ds 







we have by the change of variable formula for Riemann-Stieltjes integrals, see e.g., 
g(B t ,A t )-g(B s ,A s ) = J a 1 ' {<P{A U ,B U )) dT B u 



t rB u q, 

a"{4>(A u , v))—!-(A u , v) dv dA u , 

OX i 



s JO 



Since A' t = f(B t , A t ) we finally get 

g(B t , A t ) - g(B s ,A s ) - a"(<p(A u , v ))^-(A u ,v) dv f(B u , A u )du 

= f a'(X u )d-B u , 

J s 

which shows the assertion. 

□ 

The next Lemma can be shown by a density argument. 
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Lemma 7 Let g,h £ C([0, 1]) and denote Ah k j n = h{(k+l) / n) — h(k / n) for k = . . . , n—1, 
n £ N. // 



sup 

*e[o,i] 

as n — > oo, then it follows 



n— 1 



n 



^l [0ii] (/c/n)(A/i fc/n ) 2 -t 



fc=0 







sup 

*e[o,i] 



n-l ,. t 
i y2g(k/n)l m (k/n)(Ah k/n ) 2 - / 5 (s) ds 



as n — > oo. 



Now we finally prove Theorem Q] and [2j In the following we will denote by C random 
constants, which depend only on b, a, their derivatives, xq, H, ||-B||oo and ||-B||a with 
A < H, regardless of their value. We start with the proof of Theorem [2j 

Proof of Theorem [2j (1) We first establish a rough estimate for the pathwise error of 
the Euler scheme. For this, we follow the lines of the proof of Theorem 2.6. in |30| . Set 

A* = <P( X k/m- B k/n), k = 0,...,n 
for n E N. By a Taylor expansion, the properties of cj) and Lemma 0] we have 

Ak+l - M; = Q^^ Xk / nl ~ B k/n)( X (k+l)/n ~ X k/n) ~ ~ B k/n) AB k/n 



+ 2 ~d3^ Xk l n ' ~ Bk l n ^ a ^ Xk / n ^ 2 ^ Bk / n ^ 2 

+ n~a~~l( X k/m - B k/n)( AB k/n 



,2 



2dx 



2 



dx 2 8xi ^/"' ~ Bk /n) a ( X k/n)( AB k/n) 2 + R^ 



with 



Since 



R®\ < C({AB k/n f + A ■ AB k/n + A*). (38) 



~2 ~dxi ^ fc / n ' ~ Bk / n ^ a<T '^ Xk /n) 

= 2 Q^2 ( X k/m - B k/n)°~( X k/n) + 2 ~Q^ Xk l w ~ B k/n) ~ Q X2 Q Xl ( X tyn> - B k/n)°{ X k/n) 

by Lemma [2] (c) , we have 
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A n k+1 =A n k + b(Tl /n )-^-(xl, -B k/n )A + Q k + r£\ 



for k = 0, . . . , n — 1, n G N, with 



1 



k = 0, . . . , n. 



Qk = -2( a(J ')( X k/n)(^ B k/n) {X k , ~B k / n ) , 

Since X k / n = (f){A k , B k / n ) and using Lemma[2] (d) and (e) we get 

A n k+1 = A n k + f(B k/n ,A n k )A + Q k + R k l) 
for k = 0, . . . ,n — 1, n G N, with the function / given in Lemma [5j Note that 
sup sup \Al\ <exp(cexp(c (H^Hoo + \\B\\ 2 1/2 ))), 

neNfc=0,...,n 

as a consequence of Lemma HI Now set 

efc = A k/n -A%, k = 0,...n. 
We have eo = A — 4>(x , 0) = and 

«(fc+l)/n 



|efc+i| < \e k \(l + CA) + \Q k \ + \R i k 1) \ + 



k/n 



f(B T , A T ) - f(B k/n , A k/n ) dr 



Since 



(k+l)/n 



k/n 



f(B T ,A T )- f{B k/n ,A k/n )dr 



•(fc+l)/n 

<C I \B T -B k/n \dT + CA 2 

k/n 



we can rewrite the above recursion as 



|e fc+1 | < \e k \(l + CA) + \Q k \ + \R^\ 



with 



< C (II.BII^A 3 ^ + \\B\\ H _ £ A H+1 - £ + A 2 ) < CA H+1 - 
\Q k \ < C(AB k/n ) 2 < CA 2H ~ 2 \ 



Since also 
we get by Lemma [3] 



n— 1 n— 1 

max \e k \ < TT(1 + CA) V(|Q,-| + l^l) < CA 2 ^ 1 

:=0,...,n z — ' J 



-l-2e 



i=0 
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Moreover, we have 



max \X k/n -Xl /n \<C^ 2H - l - 2 \ 

k=0,...,n 



(41) 



due to X t = <f>(A t ,B t ), t e [0,1], and Xl /n = <f>(A%, B k/n ), k = 0,...,n. 



(2) Now we derive the exact asymptotics of the error of the Euler scheme. We can write 
the recursion for the error e k = A k / n — A k as 

e k+1 =e k + f y {B k/n , A k/n )e k A + Q k + + -f yy (B k/n , rj k )e k A 

with r\ k between A k j n and A k . Put 



Qk = -^(^^ / )(X k / n )(AB k / n ) 2 -^—(X k / n , -B k / n ). 
\Q k - Q k \ < CA^- 1 -^. 



By (HQJ) we have 
Since moreover 



(fc+l)/n 



k/n 



f y (B t ,A t ) dt - f y (B k / n , A k / n )A 



< CA H+1 - £ , 



we get by (ggj) and gU) 

r{k+X)/n 

ek+i = e k + e k / A t ) dt + Q k + R 

J k/n 



(3) 
fc 



with 

Applying Lemma [3] yields 



|^(3) | < ^,^min{4ff-l-4£,H+l-s} 



fc-1 fc-1 



i=0 j=i+l 



A k /n -A% = J2Q/l[ ( 1 + f y (B t , A t ) dt)+R 



l fc ' 



with 



sup |4 4) | 

fc=0,...,n 



sup 

k=l....,n 



fc-1 fc-1 / r fj + i)/ n 

(3) — ' 



E«f n i+ 

i=0 j=i+l \ 



j/n 



fy(B t ,A t )dt 
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Thus it remains to consider the term 

fc-l fc-i 



Now set 



Since 



i=0 j=i+l 



/ f y (B t ,A t )dt, j = 0,...,n-l. 

Jj/n 



fc-l 



'fc-l 



n ( i+a i) - ex i» n ' 



j=i+i 



< 



2-1 



exp -£\ 



3=0 



i-1 \ fc-l /fc-l 
exp | ^ a 3 - J JJ (1 + a,-) - exp 0j 

3=0 ) j=i+l \j=0 



we get 



fc-i /fc-i 
II (l + «i)-exp ^Jaj 
i=i+l \i=i 

by a straightforward application of Lemma [3l Hence we have 



sup 

0<j'<fc-l<n 



< CA 



k 1 / pk/n \ 

A k/n -A n k =Y,Qi exp / f y (B s ,A s ) ds + R™. 
i=o \ Ji / n J 



with 



sup \W k 

fc=0,...,n 



| ^(5) | < ( o^min{4i?-2-4e,f/-e} 



Moreover, since X t = 4>(A t ,B t ), t £ [0, 1] and ~X k / n = 4>(A k , B k / n ), k = 0, . . . , n we 
-X"fc/n — X k j n = -^-(A k / n , B k / n )(A k / n — A k ) + --g^(6 k ,B k / n )(A k / n — A k )' 
with between A fc/n and It follows by (00]) and (fl2|) 

-(A fc/n , fl fc/ „) ^ Qi exp / f y (B s ,A s ) ds + 
1 i=o V 1 ^/™ / 

| fl (6)| < ^j^min{4H-2-4e,H-e} 



^fc/n - ^fc/n - o 



with 
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Since finally by Lemma [6] 



dj> 
dx 



> / [ k / n \ 

-(A k / n ,B k / n )Qie-x.p / f y (B s , A s ) ds \ 

1 \Ji/n J 



"2dxi^ k ^ n ' Bk/n)(& a '){Xi/n)-Q^-{Xi/ni ~ B i/n) 



exp 



' rk/n \ 

j( f y (B s ,A s )dsj(AB i/n ) 



we have 



Xk/n ~ X k / n 



' { X i/n)Di/ n X k / n (ABi/ n ) 2 , 



. fc-l 
i=0 



Define X" = X^/ ra for t G [A;/ri, (fc + Since clearly 

we have 



|oo — max |X fc / n - X k/n \ + i? fc ' 

fe=0,...,n ' 



(7) 



with 



Thus it follows 



\r[ 7) \ < ca 



H-e 



lim n 2H ~ 1 \\X-X r 

n— >oo 



loo = lim max n 2H 1 \X k/n - Xl /n \ 

n— >oo k=0,...,n ' ' 



and we get by (|35|) 

lim n^HX-X 71 ! 

n^oo 

Furthermore it holds 



n 



2H-1 



oo = lim max — — 

n->oo fc=0,...,n 2 



fc-1 



^ v'(Xi/ n )Di/ n X k / n (ABi/ n ) 



i=Q 



ii 



2H-1 



lim max 

n— »oo k=0,...,n 2 



fe-l 



^'{X i / n )D i / n X k / n (AB i / n ) 



i=0 



lim sup 

n ^°° te[o,i] 



7i 



2#-l 



n-1 



2 l [0 , t] ^7n)(T / (X i/n )A/n^i(AS 



i/rcj 



8=0 
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with 



Z t = exp (- J b\X u ) du- J a'(X u ) cT B^j , tE [0,1]. 



This is due to the fact that the sample paths of Z are Holder continuous of any order A < H . 
It is well known that 



n-l 



n 



2H- 



k=0 



as n — > oo for all t £ [0,1]. Since n 2ff_1 ^^q 1 l[ 0)t ] (fc/n)(AB^ n ) 2 is monotone in t, the 
exceptional set of the almost sure convergence can be chosen independent of t G [0, 1]. Thus 
we get by Dini's second theorem that a.s. 



lim sup 

n ^°°te[o,i] 



n-l 



2H- 



n 



- 1 ^l [0j4] (A;/n)(Ai? fc/n ) 2 -t 



k=0 



0. 



Hence it follows by Lemma [7] 



lim sup n 2H 1 Z± 



te[o,i] 



which finally shows the assertion. 



n-l 



)D i / n X 1 (AB i/n ) 



i=0 



sup Z t 

te[o,i] 



■/ o 



a.s., 



(44) 



□ 



Proof of Theorem Q3 By ([4"5]) we have 



n-l 



Xi - = -± £ a'(X i/n )A/„*i(AB i/n ) 2 + 4 6 ) 



i=0 



with 



|.R( 6 )| < CA min{4H - 2 - 4£ ' H_e} . 
The assertion follows then by (|44|) and Lemma [71 as in the previous proof. 



□ 
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